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Abstract 

Andreev reflection at the interface between a ferromagnet and a supercon- 
ductor has become a foundation of a versatile new technique of measuring the 
spin polarization of magnetic materials. In this paper we will briefly outline a 
general theory of Andreev reflection for spin-polarized systems and arbitrary 
Fermi surface in two limiting cases of ballistic and diffusive transport. 
Andreev reflection (AR) at the interface between a superconductor and a ferromagnet 
has been attracting significant interest (e.g., Refs.0~!) as the foundation of a new technique^ 
to measure the spin polarization in ferromagnets. The technique is based on the ideal that 
the Andreev process is forbidden in half-metals, where only electrons with one spin direction 
are present at the Fermi level. Correspondingly, in a ferromagnet the Andreev current is 
partially suppressed, as not all of the conductance channels (CC) are "open" for AR; some 
channels exist in one spin direction, but not in the other, and thus do not contribute to the 
Andreev current. 

Building a stable and reliable technique for probing spin polarization based on AR is 
not always straightforward and requires a quantitative theory. Such a theory should take 
into account the following effects: (1) different number of CC for different spins, (2) finite 
interface resistance, (3) band structure effects (4) effect of an evanescent Andreev hole on 
quasiparticle current in half-metallic CC, and (5) diffusive transport in the ferromagnet, if 
needed. 
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The existing works treat only some of these questions. The first one was taken care of in 
Ref.i. The second one was answered in part in the seminal paper of Blonder, Tinkham and 
KlapwijJ, but only for "nonmagnetic" channels (the CC that exist in both spin directions). 
The third one was dealt with in Ref.i, but only in the ballistic limit. The fourth one 
was mentioned in Refs.i, but not investigated quantitatively. Finally, the last one was 
briefly touched upon in Ref.i, but left in the form that could not be directly applied to the 
experiment. Some aspects of diffusive transport and AR were addressed in Refs.S. 

In this paper, we address all these issues and present a compendium of formulas needed 
for a quantitative analysis of superconductor-ferromagnet AR. We start with a ballistic 
contact, whose size is smaller than the mean free path of electrons in the bulk. All electrons 
with a positive projection of their velocity onto the current direction, x, pass through the 
contact. Conductance of a ballistic contact is!'! 

G = j\(N\v x \)A, (1) 

where A is the contact area, N is the volume density of electronic states at the Fermi level, 
v is the Fermi velocity, and brackets denote Fermi surface averaging: 

1 „ 1 ^ r dS F 



Integration and summations are over the states with tw^ > 0, and Q is the unit cell 
volume, k, i, and a denote the quasimomentum, the band index, and the electron spin, 
respectively. It is instructive to look at Eq.|| from the "mesoscopic" perspective, using as a 
starting point the Landauer formula for the conductance of a single electron!, Go = e 2 /h. 
The total conductance is equal to Go times the number of CC, N cc , which is defined as the 
number of electrons that can pass through the contact. If the translational symmetry in 
the interface plane is not violated, then the quasimomentum in this plane, ky , is conserved, 
and N cc is given by the total area of the contact times the density of the two-dimensional 
quasimomenta. The latter is S x /(27r) 2 , where S x is the area of the projection of the bulk 
Fermi surface onto the contact plane. Thus G = ^t§^w = xl (-^l^l) ^ • 
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Let us consider now the opposite limit, when the contact size is much larger than the 
mean free path. The conductance is then given by the the bulk conductivity, which is known 
from the Bloch-Boltzmann theory: 



a 



(e 2 /h)(Nv 2 x )r, (3) 



where r is the relaxation time. The Ohm's law requires that the conductance G = aA/L, 
where L is the length of the disordered region. This can be reproduced within the "meso- 
scopic" approach!, taking into account that now each CC, that is, each separate ky state, 
has a finite probability for an electron to get through the disordered region, < T < 1, and 

G = jY. T K = jJ x ^(C)/cosh 2 (L/C), (4) 

where k = {\t,,,i,a}. T K is conveniently defined in terms of the probability distribution, 
P(C), of the localization lengths, (. The cutoff A should be of the order of the mean free 
path I] in fact, A = 2/i. Ohm's law requires that G oc 1/L, thus the behavior of P(C) at 
large ( must be const /( 2 . Normalization requires that const = XN CC . Substituting that in 
EqM we get 



G 



e 2 \N cc d( _ e 2 \N ( 



h Jx C 2 cosh 2 (L/C). hL 
e 2 AX „ r dS F 



In the constant r approximation, used in Eq.|3|, the average mean free path I = 
Eia I §<,t/ / gfu«,«, thus A K = 2v KX r. Thus 

._. e 2 A _ [2dS F 2 e 2 / 2 \ A A 

la I 1 1 

In the diffusive limit the conductance is determined by (Nv 2 ) , as it should. 

The standard theory of AR (BTK)i, places a specular barrier at the interface, and 
assumes the ballistic regime and the free electron band structure in the bulk. Let us repro- 
duce the main results of the BTK paper using, instead of their derivation, the "mesoscopic" 
approach!. Probabilities of four processes must be considered: normal reflection, AR, and 



transmission into the superconductor with or without the branch crossing!! The total cur- 
rent can be written as 



( G )ns = ji: Ts{k) = ^ £(1 + A K - fl,), (6) 

where A and S are the probabilities of the normal and Andreev reflection, respectively. 

Beenakker showedi that the "Andreev transparency" , the probability of an Andreev process, 

can be expressed in terms of the normal transparency T/v of the interface. For zero bias : 

2T 2 (k)(1 + (3 2 ) 2T 2 (l + /3 2 ) 

5 /3 2 T 2 +[l + ri v ] 2 /?2 T 2 +[2-T N ]* {) 

where Tn(k) is the normal state transparency, r 2 = 1 — Tn(k) is the corresponding normal 



state reflectance, and (3 = V/y\A 2 — V 2 \ is the coherence factor. A similar formula can be 
derived for V > A. For a specular barrier, and neglecting the possible Fermi velocity mis- 
match at the interface, Tjv(/c) = 1/[1 + Z 2 ], where Z is the BTK barrier strength parameter!. 
A simple algebra shows that Eq. [F] is equivalent to the BTK formulas. 

We will now apply this approach to the diffusive AR. A diffusive Andreev contact can be 
viewed as a contact between the normal and the superconducting leads, which in addition 
to the interface, are separated by a diffusive region. The size of the region is larger than the 
electronic mean free path!. In the zero temperature and zero bias limit, Eq.|7| reads: 

e 2 e 2 2T 2 

where now the normal state transmittance for the conductance channel k is given by the 
sequential conductor's formula: 

f- 1 - 1 = (T/ - 1) + (t- 1 - 1), (9) 

where t is the transmittance of the diffusive region, and T/v is the barrier transparency. 
Using Eq. for the distribution of t's, we find 



e 2 2 
(Gns) l = -rJ2 



h V (2/Tv - 2 + 2/t K - 1) 
e 2 \N cc r°° dy 



h L Jo [2(1 -T N )/T N + cosh y} 2 ' 
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(10) 



The last integral can be taken analytically and gives 

e 2 XN rr w cosh w — sinh w 



(11) 



h L sinh 3 w 

where cosh it; = 2(1 — Tn)/Tn. For the clean (no-barrier) interface, T/v = 1, w = «7r/2, 
and this expression reduces to Eq.[|, thus reproducing the known result0'i that the diffu- 
sive Andreev contact with no interface barrier at zero bias has the same resistance in the 
superconducting and in the normal states. 

Is it possible then to distinguish between the spin-polarization suppression of the Andreev 
current and possible diffusive transport effect using the experimentally measured conduc- 
tance? The answer to this crucial question is yes, as we demonstrate in Fig. 1: although 
it is very difficult to discern the effect of a finite Z in a ballistic contact from the effect of 
diffusive transport, it is easy to separate both of them from the conductance suppression 
due to the finite spin polarization. 

We will now derive a full set of formulas for arbitrary bias, temperature, and interface 
resistance for both ballistic aad diffusive regimes, generalizig the BTK formulaj ia order 
to be able to use them for the half- metallic CC and in the diffusive limit. These general 
formulas are summarized in Table |. 

We start with extending the BTK approach over the half-metallic CC, which, by defi- 
nition, correspond to the kn allowed in one spin direction, but not in the other. Following 
BTK, we consider an incoming plane wave and the transmitted plane wave (with and without 
branch crossing) 

(A 



Akx. 



1> 



tr 



V / 











e ikx + d 


^—ikx 


V) 









assuming, for simplicity, the same wave vector for all the states. Here u and v have the 
standard BTK meaning, u 2 = 1 - v 2 = (1 + f3)/2. Unlike BTK, though, now the reflected 
state is a combination of a plane wave and an evanescent wave: 
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ikx 
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The total current is - 4(i-/32)z(^-z) 4 +[i+(^-2Z)2][(i+/3)2+4/32z2] at eV > A, where K - 
fi/k,) and zero otherwise. 

As eV A, Ghs/Go -> , and G ffS /G -> Gjv/G = 1/(1 + Z 2 ) as V -> oo. 
We will not discuss all the aspects of the non-trivial behavior of Ghs/Gq at intermediate 
biases. Importantly, Ghs/Go generally behaves non-monotonically with V, and may have 
a maximum larger than Gn/Gq at an intermediate voltage. This maximum is due to the 
fact that, although the Andreev-reflected hole does not propagate and does not carry any 
current, the Andreev process itself is allowed at eV > A and enhances the transparency 
of the barrier. This effect does not exist, though, for Z — 0, nor for K — ► oo. In the 
formulas given in Table | we used K — > oo, to simplfy the equations, since the actual value 
of K matters in a relatively narrow region of voltages above the gap. Note that the simple 
renormalization of the normal current at eV > A, used in Ref.0, gives a rather different 
result: instead of ^ 1+/3 4 p +4Z2 it gives (x+^+tz^n^z 4 ' wn i cn diverges at eV — » A + 0. 

Now we generalize the BTK formulas beyond the ballistic hypothesis. For the nonmag- 
netic CC the calculation follows Eqs.[7| and |9|. For zero temperature and a subgap bias voltage 
eV < A(T) 

( hfc/P7*(K) + [2-f N (K)}* (13) 

and 

f^ = T N l +t- 1 -l = Z 2 + t~\ (14) 

with the distribution @ for t. After some algebra we obtain 

e 2 AiV cc r (1 + ^)dy 
{Gns)l -J—Jo /5 2 + (2Z 2 + cosh 2/ ) 2 - (15) 

Factor N cc now stands for the number of CC allowed in both spin channels, A is given by 

the average mean free path for the channels in question, and thus the total conductance is 

given by (iVf 2 ), averaged over these channels. For Z = 0, this gives 



= if ("At 7 1015 



V - A 



(16) 



which starts from the normal conductivity and logarithmically diverges at V = A. For 
arbitrary Z the conductance still can be cast into an analytical form, namely 

where 



F{s) = cosh- 1 (2Z 2 + s)/ v /(2Z 2 + s) 2 - 1 

Similarly, for > A 

\ e 2 AiV ec foo 2f3dy 
{Gns)l -J—Jo /3+(2Z 2 + coshy) 2 (17) 
e 2 r /„ 2 



<<W =-^-(iV^) iT 2/3F(/3). (18) 



At Z = this reduces to 



an 



VI \ 'It A 
interesting symmetry^. At y > A we get 



/ v e2r /»r 2\ cosh" 1 (2Z 2 + 1) 

which should be used to normalize the whole conductance curve. 

Finally, for for the "half-metallic" CC, there is no conductance at eV < A. For eV > A, 

ln x e 2 XN cc r 2f3dy 

\<^HS) L 



h L Jo (/? + l) 2 + 2(2Z 2 - 1 + coshy) 



(°hs) = e -^-(Nv 2 ) i (3F[(P + l) 2 m-l], 

where the arrow in the subscript shows that these channels are allowed only in one spin 
subband. 

It is again instructive to see how this expression behaves at V ^> A : 

v e 2r /nr 2\ cosh _1 (2Z 2 + 1) 

( ggs > = TT ( Nv )i 2Z j P + i • (21) 

which is exactly twice less than the corresponding nonmagnetic limit. 



The formulas derived in this section, and summarized in Table |, finalize our task of 
generalization of the BTK equations over the finite spin polarization in both ballistic and 
diffusive limits. The finite temperatures are taken into account straightforwardly in the 
same way as in the original BTK paper and are not discussed here. We thank E. Demler 
and I. Zutic for useful suggestions. 
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TABLES 

TABLE I. Bias dependence of the total interface current in different regimes: BNM = ballistic 
non-magnetic0; BHM = ballistic half-metallic; DNM = diffusive non-magnetic; DHM = diffusive 
half-metallic. F(s) is defined in the text 

E < A E > A 

RNM 2(l+/3 2 ) 2/3 

BHM (TT^W 

DNM ^Im [F(-i/3) - F{i(3)} 2pF(J3) 

DHM pF[(l+/3) 2 /2-l]. 
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